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ABSTRACT 

We report the detection of X-ray emission from the symbiotic star V1329 Cyg with XMM- Newton. 
The spectrum from the EPIC pn, MOSl and M0S2 instruments consists of a two-temperature plasma 
with kTi ^ 0-ll-o;o2 keV and fcTa = 0.93lo;i[4 keV. Unhke the vast majority of symbiotic stars 
detected in X-rays, the soft component of the spectrum seems to be absorbed only by interstellar 
material. The shock velocities corresponding to the observed temperatures are about 300 km s~^ and 
about 900 km s~^. We did not find either periodic or aperiodic X-ray variability, with upper limits 
on the amplitudes of such variations being 46% and 16% (rms), respectively. We also did not find any 
ultraviolet variability with an rms amplitude of more than approximately 1%. The derived velocities 
and the unabsorbed nature of the soft component of the X-ray spectrum suggest that some portion 
of the high energy emission could originate in shocks within a jet and beyond the symbiotic nebula. 
The lower velocity is consistent with the expansion velocity of the extended structure present in HST 
observations. The higher velocity could be associated with an internal shock at the base of the jet or 
with shocks in the accretion region. 

Subject headings: binaries: symbiotic - stars: individual (V1329 Cyg=HBV 475) - stars: white dwarfs 
- X-rays: stars - ISM: jets and outflows 



1. INTRODUCTION 

V1329Cygni (=H BV 475) is one of a small nu mber 
of symbiotic novae (|Miirset fc NussbaumeJ Il994i ) that 
have outbursts of 3 to 7 mag and return to their pre- 
vious brightness only slowly over decade s. The canon- 
ical model (e.g. iBaratta fc Viottil 119901) for symbiotic 
novae involves a hot white dwarf star that has a sur- 
face thermonuclear flash, with the fuel supplied by a 
red giant companion. The only recorded major out- 
burst of V1329 Cyg began in about 1965 and reached 
B magnitude < 11.5 mag in October 1966. Some 
time elapsed be t ween discovery as an eru ptive object 
(jKohoutekl 119691: iKohoutek fc BossenI |1970[ ) and recog- 
nition that it is a binary, with the object first bein g 
called a proto-planetary nebula (jCrampton et al.lll970l ). 
The photographic brightness dropped to my < 16 mag 
over the past century, with several abrupt drops of as 
much as another 2.5 mag between 1925 and 1962. These 
repeated with a period of 950-959 days and were ex- 
plained as the eclipses of the hot component by th e 
red giant fSticnon et al. | Il974t iGrvgar et aFl Il979f ). 
iSchild fc~S chmid (1997) improved the determination of 
the orbital period to 956.5 days and found from polarime- 
try an inclination of 86 ±2 degrees. From optical and UV 
emission lines, orbital parameters of the hot component 
have been determined yielding minimum masses of about 
0.71 and 2 Mq for the white dwarf and red giant, respec- 
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lively (ISc hild k Schmidl[T997| ) . These authors derived an 
Eb-v value of 0.6, correspondin g to nn = 2.3 x 10^^ cm^ 
(usin g the conversion factor of iGroenewegen fc Lamera 
[l989h . 

About 200 s ymbio tic stars are known (e.g. 
'Belczynski et alj I2000D . but jets have been de- 
tected at differe n t wav elengths only in 10 of them 
(jBrocksopp et all 120041 ) V1329 Cygn i is a mem- 
ber of this list, since iBrocksopp et al.l (|2003[ ) found 
two peaks and extended emission in HST WFPC2 
snapshot images taken in October 1999 with F502N 
and F656N filters. These were separated by about 
950 AU, a ssuming a distance of 3.4 kpc as given by 
iMiirset fcNussbaumer (.1994). After c omparison with 
HST images of ISchUdASchmidl ([1991 taken in July 
1991, they derived an expansion velocity of 260 ± 50 
km s~^, suggesting that this mass ejection was not 
associated with the nova outburst in 1965, but with 
an event in 1982. As an additional riddle, the position 
angle of the orb i tal p lane of 11 ± 2° determined by 
ISchild fc SchmidI ()1997l ) furthermore suggests that the 
mass ejection occurred along the orbital plane, instead 
of perpendicu lar to it as expected for a jet. Using the 
ephemeris of S child fc S chmid (1997) 



2444890.0 -f 956.5 x E, 



(1) 



the system was at phase 0.733 during the observations of 
ISchild fc Schmi d (1997) and at phase 0.880 during that 
of IBrocksopp et al., ()2003i ) . Whether the different ap- 
pearance of the extended [OIII] emission is due to this 
slightly different orbital phase and thus a changed ioni- 
sation structure in the stellar environment, or due to real 
expansion of the gas in side a wind or jet, i s unk nown. We 
also have to note that ISchild fc Schmid (I1997D used the 
FOC camera with filter F501N. while IBrocksopp et al.l 
(|2003f) used the WFPC2 camera with filter F502N, thus 
different filters with different properties. 
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Table 1 

Observations on May 02/03, 2009 (ObsID: 0604920301) 



Instrument 


Filter 


Duration 


UT Start 


UT Stop 


pn 


Medium 


57434 




11:00:56 


02:58:10 


MOSl 


Medium 


59021 




10:38:33 


03:02:14 


MOS2 


Medium 


59027 




10:38:33 


03:02:20 


OM 


U 


lOx - 


1500 


10:47:01 


16:15:16 




UVWl 


lOx 


1500 


16:50:37 


21:48:52 




UVM2 


lOx ^ 


1500 


21:54:13 


03:04:49 



X-ray observations provide a direct probe of the two 
most important components of jet-driving systems: the 
bow and internal shocks of the jet emitting soft X-rays 
and the central parts of the jet engine, where gas is being 
accreted to power the jet, leading to hard and /or super- 
soft X-ray emission. Cur rently RAqr (Kellogg et alj 
200l [20071) and CH Cvg (IGallowav fc So koloski 2004; 



Karovska et all l20?y7l mm are the only two jets from 



symbiotic stars that have been resolved in X-rays. 
All objects with jets, detected in other wavelength, 
when observed in X-rays, show soft components with 
kT < 2 keV (additionally to sometimes presen t swjper - 
soft emission): R Aqr j KcUogg et al.' '20 01L I2007D 
CH Cvg (.Galloway &: Sokoloski ^.2004- .Karovska et al.l 
r—^ - "loL - - - 



[20071 1201 d) MWC 560 (IStute fc Sahal I2009D . RS Oph 
(iLuna et all I2009D. AG Dra dCgnzalez-Riestra et al] 
120080 ■ Z And ([Sokoloski et al.l l2006[). The three objects 
CH Cyg, R Aqr, MWC 560 also emit h a rd co mponents 
(iMukai et al.l [20071: [Nichols et all [2007[: [Stute fc Sahal 
[20091 ). Z And s howed hard emission in one of three ob- 
servations only ([Sokoloski et al.l[2006f ). 

The paper is organized as follows: in Section [21 we 
show details of the observations and the analysis of the 
data. After that we describe the results in Section |3| We 
end with a discussion and conclusions in Sections |4| and 

B 

2. OBSERVATION AND ANALYSIS 

We observed the field of V1329 Cyg with XMM-Newton 
in 2009 (Table [1]) using the EPIC instrument operating 
in full window mode and with the medium thickness fil- 
ter. Simultaneously, we used the Optical Monitor OM. 
All the data reduction was performed using the Science 
Analysis Software (SAS) software packagcl version 8.0. 
We removed event at periods with high background levels 
from the pipeline products selecting events with pattern 
0-4 (only single and double events) for the pn and pat- 
tern 0-12 for the MOS, respectively, and applying the 
filter FLAG=0. The resulting exposure time after these 
steps is 37.6 ks. 

The source spectra and light curves were accumulated 
from a circular region of 12" radius (240 pixels. Fig. [ij 
centered on V1329 Cyg using the coordinates from SIM- 
BAD. The background spectra and light curves were ex- 
tracted from a source- free region on the same chips taken 
within an circle of 30" radius (600 pixels). Spectral redis- 
tribution matrices and ancillary response files were gen- 
erated using the SAS scripts rmf gen and arf gen, and 
spectra grouped with a minimum of 25 counts per energy 
bin were fed into the spectral fitting package XSPEC0 



vl2.5.1. For timing analysis, photon arrival times were 
converted to the solar system barycenter using the SAS 
task barycen. 

3. RESULTS 

3.1. Images 

We detected V1329 Cyg with each of the instruments 
aboard XMM-Newton. The total numbers of counts in 
the EPIC pn, MOSl, and M0S2 cameras were 329, 84, 
and 107, respectively (Fig. [1]). For comparison, the total 
numbers counts expected due to background radiation 
alone were 83, 36, and 20 for the three X-ray instru- 
ments. The average magnitudes in the optical filters U, 
UVWl and UVM2 are 12.32, 12.34 and 13.35, respec- 
tively. The X-ray images of V1329Cyg show a point 
source morphology. In order to assess the possibility of 
extended emission, we constructed radial profiles from 
the images using the SAS task eradial (a routine to 
extract a radial profile of a source in an image field and 
to fit a point spread function [PSF] to it) keeping the 
centroid fixed at the source position of V1329 Cyg and 
choosing a step size of two pixels. With this method, we 
did not find any evidence for extended X-ray emission. 

3.2. Spectra 

The spectrum of V1329 Cyg is relatively soft with pho- 
tons having energies of < 2 keV (Fig. [2|). It is well 
described (x^ = 1-10, 6 d.o.f.) by a two-temperature 
plasma model with N overabundant by a factor of ~ 
4 in the softer compone nt compared to solar values 
([Anders fc Grevessd [198 91) in order to explain the high 
flux around 0.4 keV (vapec in XSPEC). For the soft 
component, the flt is consistent with no internal ab- 
sorption, with Nh in the range (0 — 5) x 10^^ cm~^. 
For the hard component, we find Nh = 9^0 32 x 10^^ 
cm^^. These values include internal absorption as well 
as interstellar absorption which in the direction of the 
source is about 2.3 x 10^^ cm~^. The fit is slightly 
worse (x^ = 1.27, 4 d.o.f.) including (non-vanishing) 
absorption. The temperature of the cooler component is 



kTi = 0.11 



-1-0.02 
-0.02 



keV, while the temperature of the hot- 
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ter component is fcTz = 0.93to;i4 keV. The fiux of the 
hotter component is three times lower than that of the 
soft component. The total model fiux between 0.3 and 
2 keV is 1.21 x 10~^* erg s~^cm~^. With the assumed 
distance of 3.4 kpc, this fiux corresponds to a luminosity 
of 1.7 X 10^^ erg s"^ 

The slope of the ultraviolet spectral energy distribution 
indicates that the optical fiux is dominated by nebular 
emission and the red giant (Fig. [3]). The flux level of 
blackbody emission corresponding to a white dwarf with 
a radius of 8 x 10^ cm at a dista nce of 3.4 kpc with an 
effective temperature of 145000 K ([Miirset et al.l[1991[ l is 
lower than the measured fluxes by 1-2 orders of magni- 
tude. 

3.3. Light curves 

We examined time series binned at 3000s in the energy 
range 0.3-2 keV (Fig. Hj) and the OM photometry with 
exposure times of about 1500s (Fig. [5]). We chose the bin 
size for the X-ray light curve to produce at least roughly 
25 counts per bin while still providing time resolution 



Figure 1. Left: XMM-Newton EPIC pn image (detail) of tiie region around V1329 Cyg in the 0.3-2 keV range. Also shown is the 
extraction region for the source events. Contour levels are 2, 3, 4 and 5 counts per pixel; middle: XMM-Newton EPIC MOSl image of the 
same region. Contour levels are 1.5, 2 and 2.5 counts per pixel; right: XMM-Newton EPIC MOS2 image of the same region. Contours are 
1, 2, 3 and 4 counts per pixel. 

mean value. Only in filter UVWl, the ratio s/sexp > 1 
indicates variability. 

We also searched for modulated emission in the light 
curves. Using the arrival times of the source events, 
we ca lculated the Zf (Rayleigh) statistic (jBuccheri et al.l 
119830 in the frequency range fmin = l/TLxp to /max = 
l/2tframe with A / = l/Texp, where Texp is the effective 
exposure time of 37.6 ks and tframe is the readout time 
(73 ms in the case of XMM). The value of needed to 
detect pulsation with a probability P is given by 



10-» 



Energy (keV) 

Figure 2. Observed spectrum of V1329 Cyg together with a model 
consisting of two absorbed thermal components (vapec; dotted 
lines) with variable abundances. 
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Figure 3. Spectral energy distribution in the optical filters U, 
UVWl and UVM2. Also plotted is blackbody emission represent- 
ing a white dwarf with a temperature of 145000 K and radius of 
8 X 10* cm. The slope and flux level of the spectral energy dis- 
tribution indicate that the optical flux is dominated by nebular 
emission and the red giant. 

of less than an hour. The X-ray and two of the optical 
light curves are consistent with a constant flux. The 
measured rms variations, those expected from Poisson 
statistics, and the ratios of these two quantities are listed 
in Table [2j s and Sexp are given in percentage of the 



Z^>2x ln(. 



Toxp X A/ 



P 



(2) 



Thus for a 3-cr (P = 2.699 x lO"^) detection, Z^ > 36.7 is 
required. The highest value of Zf detected in our obser- 
vation is 27.5, which corresponds only to 1.26-a. With 
A^s=329 and A^_b=83 source and background photons de- 
tected with the EPIC pn camera respectively, the pulse 
fraction corresponding to the highest Z^ value detected 
is 



p = PobsiNs+NB)N. 
^[2Zf{Ns + NB)V^^N. 
= 0.46 ±0.009 



±[2{Ns + Nb)]^^^N. 



-1 
s 



(3) 



for a nearly sinusoidal signal, where po^s is the pulsed 
fraction uncorrected for the background contribution. 
Therefore, our observation is only sensitive to pulse frac- 
tions p > 46 %. 

4. DISCUSSION 

The X-ray images of V1329 Cyg show a point source 
morphology; no hints of extended emi ssion are present. 
In Fig ure El we show HST images from iSchild fc Schmidl 
(fT997[) and lBrocksopp et al.l (|200l where extended emis- 
sion in [OIII] is evident. These authors concluded that 
the visible optical emission is an expanding jet, however, 
resulting in the conundrum that the mass ejection oc- 
curred along the orbital plane. 

The two-temperature X-ray emission can be exp lained 
in the light of simulat i ons o f jets in symbiotic stars (jStutd 
120061 iStute fc Sahail 120071 ) where soft X-ray emission 
arises from internal shocks and the bow shock. The ve- 
locity of the shock, Ushock, can be derived using the mea- 
sured temperature of the component and the following 
relation (assuming strong shock conditions): 



post shock — w ^ 



3 /i mp 

16 



■ V, 



sliock 
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Table 2 

Measured and expected variations and their ratio 



Instrument 


measured variation s 


expected variation Sexp 


ratio s/sexp 


X-ray 


39.62 % 


38.56 % 


1.03 


optical U 


0.054 % 


0.061 % 


0.89 


optical UVWl 


0.095 % 


0.025 % 


3.77 


optical UVM2 


0.063 % 


0.067 % 


0.94 



0.010 - 




0.002 



0.000 



4x10' 



Figure 4. EPIC pn light curve binned at 3000 s in the energy 
range 0.3-2 keV. The dashed line shows the average count rate. 
Overall the light curve is consistent with X-ray flux being con- 
stant - the ratio of measured rms to that expected from Poisson 
statistics, s/sexp, is 1.03. 



: 0.105 keV 



'^shock 



.300 km s- 



(4) 



with nip the proton mass, ks the Boltzmann constant 
and fi = 0.6 the mean particle weight. Therefore the soft 
component temperature corresponds to a shock velocity 
of about 300 km s~^, consistent with the expansion ve- 
locity of 260i b50 km s~^ of the optica l extended structure 
measured bv iBrocksopp et all (|2003[ ). 

Furthermore the temperature of the hot component 
suggests the presence of another shock with a velocity of 
about 900 km s~^. The spatial resolution is by far not 
good enough to further connect the apparent flow in HST 
observations with our results. However, if the optical ex- 
tended structure represents a jet, then the second shock 
could be located at the base of the jet. The harder of the 
two X-ray components could also arise from the bound- 
ary layer between th e accretion disk and the white dwarf. 
iPandel et al.l ()2005D found that the maximum temper- 
ature in a cooling flow representing an optically-thin 
boundary layer is fcTmax = 3/5 fcTviriai- The tempera- 
ture inferred from a single-temperature flt will be lower 
than fcTinax- The plasma temperature fro m single tem- 
peratu re fits to boundary layer emission in lKennea et al.l 
(l2009l) from RT Cru, T CrB, CH Cyg and SS73 17 are 
lower than the maximum te mperatur es determined from 
cooli ng flow fits in RT Cru dLunafe Sokoloski 2007), T 
CrB (ILuna et al.l[2008h and SS73 17 ( Eze et al.li201Q) . 
The higher absorbing column in front of the hard com- 
ponent is consistent with this component being closer to 
the central engine than the lower temperature plasma. 

The fact that the soft component of the X-ray spec- 
trum is consistent with no internal absorption can only 
be explained if the source of the X-ray emission lies out- 
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Figure 5. Light curves of the OM exposures for filters U (top, 
3440 A), UVWl (middle, 2910 A) and UVM2 (bottom, 2310 A). 
The dotted lines show the average magnitude. The light curves for 
filters U and UVM2 are consistent with a constant flux {s/sexp = 
0.89 and 0.94), while s/s^xp = 3.77 for the UVWl fllter indicates 
variability. 

side of the symbiotic nebula with a typical size of a few 
AU. In contrast, almost all symbiotic stars detected in 
X-rays show spectra absorbed b y columns densities of a 
few time 10^^ to 10 ^^ cm^ (e.g. iLuna fc Sokoloskil [20071 : 
iKennea et ani2009f ). 

The over-abundance of N by a factor of « 4 has 
some implications for the o rigin of the shocked mate- 
rial. ISchmid fc SchiidI (Il99flh find in V1329 Cyg an en- 
hancement factor of about 10 in nitrogen in the symbi- 
otic nebulosity from optical, NIR and UV spectroscopy. 
The mean abundances found in M giants show an over- 
abundance of N by a factor of 3.5 compared to solar val- 
ues (jSmith & Lambert 1985, 1986), due to the conversion 
of carbon into nit rogen and mixing into th e atmosphere 
of the red giant (|Nussbaumer et ahl 119881 ). Larger val- 
ues may be attributed to more advanced nuclear burning 
stages. Therefore the shocked material in the softer com- 
ponent was probably initially part of the red giant wind. 
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Figure 6. Left: HST image in tlie filter F501N of tfic region around V1329Cyg from ISchild fc SchmidI l(T997l). The orientation of the 
orbital plane is indicated by the white line; right: HST image in the filter F502N of the same region from IBrocksopp et al.l pOOl ). For 
comparison, a scale of 300 mas has been added in both images. 



which has not been further reprocessed during thermonu- 
clear burning close to the white dwarf (jSchmid &: Schildl 
[1990) and has been ejected outside of the symbiotic neb- 
ula within a jet. 

The possible lack of short term variability in the X- 
ray light curve may support the above scenario. Vari- 
ability with time scales of minutes to hours is typically 
observed in symbiotics where X-ray originate in the ac- 
cretion flow (e.g. SS73 17, Eze et al. 2010, or CH Cyg, 
iMukai et al.ll2007|) . However, because the X-ray count 
rate from V1329 Cyg leads to large error bars, we cannot 
definitely exclude flickering. Variability is statistically 
detected only in the UVWl band at 2910 A and the rms 
amplitude of this variability is only 1%. 

5. CONCLUSION 

We report the detection for the flrst time of X-ray 
emission from the symbiotic star V1329 Cyg which was 
not previously detected by ROSAT due to its low X-ray 
flux. The observed X-ray temperatures and especially 
the unabsorbed nature of the soft component of the X- 
ray spectrum indicate that some of the X-ray emission 
might originate in shocks inside a jet outside the symbi- 
otic nebula. In this regard, VI 329 Cyg adds to the small 
but growing group of jet-driving symbiotics with X-ray 
emission. 
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discussions. This work is based on observations obtained 
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